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Equations of Maxwell Type 

K. O. Makhmudov, O. I. Makhmudov, and N. Tarkhanov 



Abstract. For an elliptic complex of first order diff'erential operators on a 
smooth manifold X, we define a system of two equations which can be thought 
of as abstract Maxwell equations. The formal theory of this system proves to 
be very similar to that of classical Maxwell's equations. The paper focuses 
on boundary value problems for the abstract Maxwell equations, especially on 
the Cauchy problem. 
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Introduction 

The term Maxwell's equations applies to a set of eight equations published by 
Maxwell in [Max65j . He called them "general equations of the electromagnetic 
field." 

To be more specific, we will restrict ourselves to reduced Maxwell's equations 
which describe the propagation of electromagnetic waves in a homogeneous isotropic 
medium in M.^ with electric permittivity e, magnetic permeability /i and electrical 
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conductivity a. An electromagnetic wave is described by the vectors of electromo- 
tive force E{t,x) and magnetic force H{t,x). Then, Maxwell's equations are 

-'e{d/dt)E + curl H = aE, , 
n{d/dt)H + c\\y\E = 0, ^ ' 

see for instance |Jac75| . 

Maxwell's equations have a close relation to special relativity. Not only were 
Maxwell's equations a crucial part of the historical development of special relativity, 
but also special relativity has motivated a compact mathematical formulation of 
Maxwell's equations. 

Although Maxwell's equations apply throughout space and time, practical prob- 
lems are finite and solutions to Maxwell's equations inside the solution region X are 
joined to the remainder of the universe through boundary conditions, cf. MonOSj, 
[SV95| . |Fri90j . etc., and started in time using initial conditions, cf. for instance 
|HH94| . 

In [MNT08] we constructed an explicit formula which restores solutions of the 
oscillation equation of the couple-stress theory of elasticity in a bounded domain X 
in through the given displacement and stress values on a part S of the boundary 
of the domain. The main ingredient of our construction is an expansion of the 
fundamental solution <P{x — y) of couple-stress elasticity. The equations of couple- 
stress elasticity factorise the Helmholtz equation in M.^ . It follows that — y) 
amounts to the quotient applied to the fundamental solution of the Helmholtz 
equation. In order to expand this latter we used spherical harmonics in R"^ and 
Bessel functions. 

Maxwell's equations factorise the Helmholtz equation in R'^, too, and so the 
construction of [MNTOS] goes through in this case in much the same way. This 
work was intended as an attempt at constructing an explicit formula which restores 
solutions of Maxwell's equations in a bounded domain X in 'M? through their values 
on a part S of the boundary of X . While such a formula is of practical interest 
and, to our best knowledge, new, it would not be surprising, for the arguments are 
routine. 

However, on thinking of the so-called Stratton-Chu formulas jStr41j we have 
revealed a very rich algebraic structure which goes beyond classical Maxwell's equa- 
tions. Eliminating the electric or magnetic field from Maxwell's equations leads to 
the Helmholtz equation for the components of the other field. Algebraically the 
elimination means multiplication of Maxwell's operator with a matrix of scalar dif- 
ferential operators from the left. Any such matrix is referred to as the quotient 
of the Helmholtz and Maxwell's operators. In contrast to Dirac type operators 
which multiply by their formal adjoints to the Laplace operator, the quotient of the 
Helmholtz and Maxwell's operators is of order 2. 

Hence it follows that the potential theory for Maxwell's equations is not quite 
standard. To handle simple and double layer potentials related to these equa- 
tions one has to take care of appropriate modification of the results in |Tar90| 
Ch. 2]. The nature of such modification becomes most transparent if one specifies 
Maxwell's equations within a wider class of differential equations associated with 
elliptic complexes. 

For this purpose, we recall in Section [1] the concept of Maxwell's equations. 
When written in terms of differential forms, they can readily be interpreted in the 
context of arbitrary elliptic complexes. The equations obtained in this way we 
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call abstract Maxwell equations, and they raise new boundary value problems. In 
Section [2] we generalise the Stratton-Cliu formula to solutions of abstract Maxwell 
equations. This allows one to construct explicitly the so-called Calderon projec- 
tions, see Theorem 5.4 of ^PS02i 1.4.1], which reduce the study of boundary value 
problems to boundary integral equations, cf. Section [3l The Cauchy problem with 
data on the whole boundary for classical Maxwell's equations is certainly overde- 
termined. In Kra 92] and [KKR03| . some solvability criteria for this problem are 
proved by using biquaternionic functions. The Cauchy problem with data on a 
nonempty part S of the boundary is of great importance for analysts. We discuss 
this problem in Sections H [5] and IB The method of [YIM92| and [Mak04| apphes 
to construct a Carleman function for the Helmholtz operator, and so for all its ma- 
trix factorisations, in the case of conical domains. Here, we use another approach 
which is elaborated in [Tar95j and Shl96, . In Section [7] we apply these results to 
classical Maxwell's equations. 

1. Abstract Maxwell's equations 

For the sake of simplicity we confine ourselves to the case of time-harmonic 
electromagnetic waves. Then the electric and magnetic fields are of the asymptotic 
form 

E{t,x) = ie + ia/u;)-^/^e-"^*E{x), 
H{t,x) = fi-^/^e~"^*H{x). 

From Maxwell's equations (lO.ip with arbitrary time dependence we deduce easily 
that the parts of E and H depending on the space coordinates satisfy Maxwell's 
equations 

tkE + cmlH = 0, 
-^fc H + curl E = 0, 

where k is the wave constant defined by /c^ ~ (e + ia/LL>)^uj'^. The sign of k is chosen 
from the condition 3fc > 0. 

The expression of Maxwell's equations in terms of differential forms leads to 
a further notational and conceptual simplification. On introducing the de Rham 
complex in R'^ 

^ /2i(R3) n^(R^) 

we can think of i? as a differential form u of degree 1 , 7J as a differential form / of 
degree 2, thus identifying curlE' with du and mrlH with d*f. Here, d* stands for 
the formal adjoint operator of d. In this way Maxwell's equations can be written 
in the form 

tku + d*f = 0, . 
-ikf + du = 0, ^ 

which already make sense not only for differential forms u and / of degree 1 and 2 
in R'^, respectively, but also for differential forms u and / of degree i and i + 1 in 
M", where —l<i<n. 

For a recent account of mathematical treatment of the electromagnetic scatter- 
ing theory we refer the reader to |CK98) . [PS02j . etc. 

Equations (jl.ip generalise to arbitrary complexes of differential operators on a 
C°° manifold X with or without boundary. More precisely, consider a complex of 
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first order differential operators on X acting in sections of vector bundles over X, 
i.e. 

— > C°"{X,F") C°°{X,F^) ^ . .. C°^{X,F^) — > (1.2) 

where A' e BiS\F\ F''+^) satisfy A'+^A' = for all i. As usual, we write A'u 
simply Au for u G C°°{X, F*), when no confusion can arise. 

By F^ is meant a smooth vector bundle of rank fc; which is nonzero only for 
i = 0, 1, . . . , TV. We give F^ a Hermitean metric, i.e. scalar products x i— s- (v, w)x in 
the fibres F^ which smoothly depend on the point x ^ X. This defines a conjugate 
linear isomorphism * of F' to the algebraic dual bundle F" by {*w,v)x = {v,'w)x 
for v,w £ F^. 

Fix a smooth positive volume form dx on X. This yields a scalar product on 

Cc^mp('^,^^)by 

/ {u{x),v{x))xdx 



JX 

for M, u e C^j,jp(A:', F*). The completion of this space with respect to the corre- 
sponding norm is denoted by L^{X, F'). We moreover introduce the formal adjoint 
j4" G Diff (F'+^,F*) for each operator A^ by requiring {Au^g) = {u,A*g) for all 
u G C°°(A', F') and g G C°°{X.F^^^) whose supports do not meet each other on 
the boundary of X . 

The formal selfadjoint operators Z\* = A^* A^ + A'-^^A^^^* on X are called the 
Laplacians of (|1.2p . The ellipticity of this complex joust amounts to the fact that 
each Laplacian Z\* G Diff2(F') is an elliptic operator of order two, cf. |Tar90[ 
2.1.4]. 

Given any — 1 < « < iV, by Maxwell's equations for complex (|1.2p at step i are 
meant 

—e{d/dt)u + A*f = au, „^ 
fi{d/dt)f + Au = 0, 

with u and / being unknown functions of t with values in sections of F' and F*''"^, 
respectively. When looking for time-harmonic solutions to these equations of the 
form 

u{t,x) ^ {e + ta/Lu)~^^^e~"^*u{x), 
f{t,x) = A^-i/2e— 

one arrives at stationary equations for the parts of u and / that depend only on 
the points of the underlying manifold X . These are 

iku + A*f = 0, 
-ik f + Au = 0, 

where k is the wave constant defined above. 



Definition 1.1. Let — 1 < i < A^. By the Maxwell operator for complex (|1.2p 
at step i is meant 

ik A* 



/ ik A'" \ 
\ A -ik )' 



As is usual in homological algebra, we will omit the index i of AF when it is 
clear from the context. 

By definition, ilF is a first order differential operator from sections of F^QF'"*"^ 
to sections of the same bundle over X. This operator fails to be elliptic of order 1 
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in the classical sense unless N — 2. On the other hand, applying A* to both sides 
of iku + A* f = we conclude that A^^^*u = unless k — Q. Analogously, from 
—ik f + Au = it follows that A'^'^^f = unless fc = 0. Complementing Maxwell's 
equations by their differential consequences A^~^*u = and = yields a 

system of first order differential equations for u and /, whose classical symbol is 
injective. Another way of stating this is to say that there is a differential operator 
from sections of F* © F^+^ to sections of the same bundle, such that C*M* is 
a second order differential operator on X elliptic in the classical case. An easy 
computation shows that 

ik + {l/ik)A'-^A'~^* A" \ , . 

A' -ik-{l/ik)A'+^*A'+^ )■ ^ ' 



Lemma 1.2. As defined above, satisfies 

e 



Proof. This is straightforward. □ 
2. Stratton-Chu formula 

For classical Maxwell's equations, the so-called Stratton-Chu formula |Str41] 
lies in the base of application of analytical methods, cf. |CK98| 4.2]. From the 
point of view of contemporary analysis this is just a very particular case of Green 
formulas, cf. [Tar90[ 2.5.4]. 

Let A" be a compact manifold with boundary smoothly embedded into a larger 
C°° manifold A", and fO]) be defined on aU of A". Suppose that both - fc^ 
and — k^ are elliptic of order two and satisfy the uniqueness condition for the 
local Cauchy problem on A". Then these operators have left fundamental solutions 
on A" which we denote by G* and G'''^^, respectively. These latter are classical 
pseudodifferential operators of order —2 acting on sections of vector bundles 
and over A". 



Lemma 2.1. The pseudodifferential operator 

G' {tk + {l/ik)A^-^A^~^*) G^ A" 

G*+i A' G*+i {~ik - {l/ik)A'+^*A'+^) 

is a left fundamental solution of the Maxwell operator N'P on A"'. 

Proof. From Lemma [TT2] it follows immediately that 

is a left fundamental solution of A'P . It remains to substitute the explicit expression 
(O) for C\ □ 



Note that is a pseudodifferential operator of order (not —1) on sections of 
vector bundle ® over X' . This corresponds to the fact that the Maxwell 
operator is elliptic in the sense of Douglis-Nirenberg with weights Si = tj = 
provided /c ^ 0. As usual, we use the same letter to designate the operator acting 
from smooth sections of compact support of the bundle F"^ © to distribution 
sections of the same bundle, and its Schwartz kernel on A" x X' which is denoted 
by $'{x,y). 



6 



K. O. MAKHMUDOV, O. I. MAKHMUDOV, AND N. TARKHANOV 



Theorem 2.2. For each u G H'^{X,F') and f £ H^{X,F'+^) satisfying 
M^{u,f) = in the interior of X , it follows that 

^ '^'\a{A-){y,iv{y))u{y) ) \ ^ ^ ^ y, ^ y 



0, if xeX'\ X. 

Here, v{y) is the outward unit normal vector to the boundary at a point y G dX ^ 
ds the surface measure on X, and cr(^*) the (classical) symbol of . 

Proof. This formula is a particular case of the general Green formula |Tar90[ 
2.5.4], for 



□ 



Write t{u) and n{f) for the tangential part of u and the normal part of / on 
the boundary of X relative to complex (|1.2p . respectively, see |Tar901 3.2.2]. Then 
we obtain 

^a{Ar{-Mf\ ^ MA%,i.)n{f)\ 
a{A'){-,ii^)u J \ia{A^){-,i^)t{u) J ^'^ 

on dX. A familiar argument shows that Theorem 12.21 still holds for all solutions 
{u, /) of Maxwell's equations in the interior of X, such that t{u) and n{f) have 
weak limit values on the boundary dX, cf. [Tar95[ 9.4]. 

For classical Maxwell's equations the formula of Theorem 12.21 is known as the 
Stratton-Chu formula, see [Str41| . In that case it manifests a more refined struc- 
ture, for the fundamental solution can be written explicitly. This can indeed be 
done in the context of the so-called complexes of Dirac type in K" which are char- 
acterised by the property that all the Laplacians are diagonal operators, with 
the usual (nonnegative) scalar Laplace operator on the diagonal. In other words, 
Z\* = —Ek^A holds for all i from through iV, where Ek^ is the unity {ki x ki)- 
matrix and A = + . . . + d^. Let e{x) be the standard fundamental solution of 
convolution type for the Helmholtz operator A + k'^ in R", that vanishes at oo. 
E.g., if n = 3 then 

, , —1 expfifcja;]) 
^ ' 47r la;] 

away from the origin. 

Corollary 2.3. Let il.2\) be a complex of differential operators with constant 
coefficients of Dirac type, and X a closed bounded domain with smooth boundary in 
R". Ifu G H^{X,F') and f G H^{X,F'+^) satisfy M\uJ) = m the interior of 
X, then 

il/ik)A"A^ -A" ^f ( inif) 

-A^ ~{^/^k)A'A" ) 7,^'^^ \^aiA%u)t{u))''' " [fix]) ^"''^ 

for all X £ X \ dX , and the left-hand side vanishes away from X . 



Formula (|2.2p is in terms of the only operator A^ which enters into the defi- 
nition of Maxwell's equations. It extends modulo smoothing operators to general 
Maxwell's equations. 
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Proof. For a Dirac type complex one can choose G^u — ~e*u for distributfon 
sections u of with compact support in R", where by e*u is meant the convolution 
of e and u. The Schwartz kernel of is obviously Ek.e{x — y) . Hence, the Schwartz 
kernel of is 

~G'\y, dy) - y)) = ~G\y, d,)e{x - y), 

the variable y of C* can be neglected, for the coefficients are constant. From ()1.4p 
one easily deduces that 

V A' ^{l/ik){A'+^ -A'A'*) 

and so 

-(l/tk)A"A^e{- ~y) A"e{-~y) 

A^e{--y) {l/tk)A^A"el ~y) 

on M"\{y}, for ( A + k'^)e = away from the origin. To complete the proof it suffices 
to use Theorem 12.21 and observe that a{A'^){i/) = \v\^E}:. = Ei^. in our particular 
case. □ 



As is easy to check 
M' 



I 1 _^ 

ik 

-A' — -A'A' 



I -^A"(Z\*+i -fc^) 

^ ik 



\ — -^^ / V —AUA'-k^) 
\ ik / \ ik 

hence the left-hand side of (j2.2p satisfies Maxwell's equations in R" \ dX for all 
integrable functions t{u) and n{f) on the boundary. In this way we arrive at what 
is usually referred to as the Cauchy type integral associated to Maxwell's equations, 
cf. [TarQOl 3.2.3]. 

3. Boundary value problems 



The Stratton-Chu formula of Theorem 12 . 21 manifests the Cauchy data on dX of 
sections u e H^{X,F^) and / G _ff with respect to Maxwell's operator 

M*. These are t{u) and n{f). 

Lemma 3.1. Suppose fc^ is not real. Let u G C'^{X,F''-) and / e C^(A', 
satisfy M^{u, f) — in X . Either of the conditions t{u) = and n{f) = on dX 
implies {u, f) = in all of X . 

Proof. Suppose t{u) = on dX. Then {Au,Au) = {A*Au,u). The second 
Maxwell equation yields A*Au — {ik)A*f in X. By the first Maxwell equation, 
{ik)A* f = k^u in X. Hence [Au^Au) = k'^{u,u). Since the imaginary part of fc^ 
does not vanish, we get u = in X. It follows that / = in ^, showing the first 
part of the lemma. The proof of the second part is similar because of the obvious 
symmetry. □ 

Both t{u) and n{f) are sections of a vector bundle F/ over the boundary. This 
is the kernel of the bundle homomorphism a{A'^~^)* {■, v) acting from F'^\gx to 
F'-^\qx, cf. |Tar90l 3.2.2]. 

It is not our purpose to study general boundary value problems for solutions 
of Maxwell's equations. We merely discuss a particular problem which stems from 
scattering of incident electromagnetic waves by a perfectly conducting body. In 
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this case the tangential component t{u) of the "electric" field u must vanish on the 
body surface dX. 

Given a section uq of over dX, we consider the problem of finding a solution 
{u, /) of Maxwell's equations in the interior of X, such that the tangential part 
of u on the boundary is well defined and t{u) — uq on dX. If a solution exists 
and is sufficiently smooth in X, then A*u = and so = A*Au = k'^u in 
X. The Dirichlet problem for solutions of the elliptic equation (Z\' — k^)u — 
in X requires not only given values t{u) on the boundary dX but also n{u). The 
homogeneous formal adjoint problem fails to have nonzero solutions unless (and 
so fc) belongs to a discrete set of nonnegative numbers with the only accumulation 
point at infinity. In this latter case the Dirichlet problem is solvable if and only if 
the pair (uo, n{u)) satisfies some orthogonality conditions. In any case the canonical 
solution to the Dirichlet problem is written in the form u — p^{uo,n{u)), where 
is the so-called Poisson operator. On having granted the component u in A" we 
define the component / from the second Maxwell equation by / = {l/ik)Au. In 
order that the first Maxwell equation be satisfied it is necessary and sufficient that 
A*u ^ 0, for then 



in X. Summarising we conclude that the above boundary value problem is solvable 
if and only if A*p^{uo,n{u)) = holds in X. Since A and A* commute and so 
- k'^)A*u = in the interior of X, the equality A*p\uo, n{u)) = is fulfilled 
in X if and only if it is fulfilled on the boundary. In this way we arrive at the 
following result. 

Theorem 3.2. As formulated above, the boundary value problem is equivalent 
to the equation A*p^(uQ,n(u)) — on dX. 

For a treatment of this boundary integral equation in electromagnetic scattering 
theory we refer to [PS02] . 

To demonstrate Theorem 13.21 we look more closely at the extreme cases i = 
and i — N. For i = 0, the operator A'^^ is zero, and so the normal part n{u) 
vanishes and the equation A* p^{uo,n{u)) = on the boundary is automatically 
fulfilled. The theorem says that the boundary value problem in question has a 
solution for each uq G 'D'{dX,F^), and the canonical solution is actually given by 
u = p'^{uo). If i = N, then the operator A^~^* has injective symbol whence the 
bundle is zero. Thus, the boundary condition t(u) — uq is empty. However, the 
section is zero and so / = 0, which implies u = in X. This corresponds to 

the fact that the boundary integral equation A*p^{u) = for a section u of F^\gx 
has the only solution u = 0. 



When combined with Maxwell's equations, the data t{u) and n{f) actually 
determine both u and / on the boundary, for 



tku + A*f 



iku+ [l/ik)A*Au 
{l/ik){A' - e)u 




4. The Cauchy problem 




-{l/ik)n{A*f), 
{l/ik)t{Au) 



EQUATIONS OF MAXWELL TYPE 



9 



while n{A* f) and t{Au) are uniquely determined by n{f) and t{u) (the so-called 
tangential operators on dX, cf. |Tar90[ 3.1.5]). Hence, the Cauchy problem for 
Maxwell's equations in X consists in finding a solution (u, /) to M^{u, /) = with 
prescribed data t{u) and n{u) on a part of the surface dX . For a general theory 
the reader can consult |Tar95j . To obtain constructive results we restrict ourselves 
to Dirac type complexes of differential operators with constant coefficients in M", 
as defined in Section [2l 

Let S be an open piece on the boundary of X . The Cauchy problem for solutions 
of Maxwell's equations in X with data on S consists in the following. Given sections 
uo and /o of the bundle Fl over 5, find a solution (w, /) to A'P{u,f ) = in the 
interior of X, such that t{u) ~ uq and n{u) — fo on S. To study this problem 
we have to choose function spaces for uq, /o and {u,f). It is a little cumbersome, 
for the behaviour of the solution near the boundary of S requires a careful study. 
In order to highlight principal difhculties in the Cauchy problem we restrict our 
attention to the case where uq and fo are continuous sections of over the closure 
of S. If exists, the solution (w, /) should be continuous up to S in the interior of X, 
and we avoid discussion of weak boundary values of (u, /) on S. We thus consider 
the problem 

M\uJ) = in X\dX, 

t{u) = Mo on iS, (4.1) 
"(/) = /o on S. 
It is well known that this problem has at most one solution in any reasonable 
space of functions on X . 

To study ()4.1|) we introduce an integral completely determined by the Cauchy 
data uq and /o on 5, namely 

G{uo,h){x)= I #(^--)f . ]ds (4.2) 

for X G K" \ iS. Since the fundamental solution #*(x — y) is real analytic away from 
the diagonal of K" x R", it follows that G{uo, /o) is real analytic in the complement 
of 5. 

Moreover, the potential ()4.2p satisfies Maxwell's equations M''Q{uo, fo) = in 
\ S. In particular, the components of the vector-valued function G{uo, fo) prove 
to be solutions to the (scalar) Helmholtz equation {A + k'^)Q = in the complement 
of S. 

When X crosses the hypersurface 5, the integral Giuo, fo){x) has jump. The 
corresponding jump formulas look very like the classical Sokhotskii-Plemelj formu- 
las, cf. [Tar90[ 3.2.3]. 

Theorem 4.1. In order that there be a solution [u, f) of the Cauchy problem 
{4-1^ continuous up to S, it is necessary and sufficient that the integral G{uo,fo) 
might be extended from R" \ X through S to the interior of X as a real analytic 
function. 

Proof. Necessity. Suppose there is a solution {u, f) of the Cauchy problem 
(|4.ip which is continuous up to S. Define sections U and F of F^ and in 
R" \ dX by 

U\ ^ I e(Mo,/o)- (^^) in X\dX, 
1 G{uoJo) in W'\X. 
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Write for the restriction of U to the open sets X \ dX and M" \ X, respectively, 
and similarly for F . 

We can assume, by shrinking X if necessary, that both u and / are continuous 
up to the boundary of X. Using Theorem 12 .21 we get 



for all X in the interior of X. Hence it follows that C/+ and extend through 
S to real analytic sections of F' and on all of (R" \ dX) U S. Furthermore, 



applying Theorem 12.21 once again yields 



dx\s 



' ia{A')[v)t{u) J Js \ia{A'){v)t(u) 



for a; £ M" \ X. Therefore, Q{uq, fo) extends from M" \ X through S to the interior 
of A" as a real analytic function, as desired. 

Sufficiency. Conversely, let U and F be real analytic sections of bundles 
and over (R" \ dX) U 5, such that 

^ =G{uoJo) 

away from X. Then AP{U,F) = in M" \ X. Since ]VP{U,F) is real analytic in 
(M" \ dX) U 5, it actually vanishes in X, too. 
Set 

u{x)\ . . (U{x) 

for X in the interior of X . From what has already been proved it follows that u and 
/ are continuous up to S and satisfy Maxwell's equations. We claim that (m, /) is 
the desired solution of (|4.ip . To see this, it remains to verify that t{u) — uq and 
"-(/) = fo on S. 

To this end, denote by Geiuo, fo) and Qmiuo, fo) the components of Q{uo,fa) 
that are sections of the bundles and F'+^, respectively (electric and magnetic 
fields). Since U and F are continuous in (R" \ dX) U S, we get by the Sokhotskii- 
Plemelj formulas 

t{u) = t{g,{uo,fo)+ -u+) 

= t{ge{uoJo)+ -U-) 

= t{ge{uo,fo)+~ge{uo,fQ)~) 
= Uo 

on S. Analogously, 

n{f) = n{grn{uoJo)+ -F+) 

= n (Gmiuo, fo)^ - F^) 

= n(f/„(wo,/o)^ -5m("o,/o)") 

- fo 

on S, which is our assertion. □ 
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Theorem 14.11 may be summarised by saying that any expUcit formula for ana- 
lytic continuation from R" \ X through S to the interior of X leads to a formula 
for solutions of the Cauchy problem with data on S. This idea goes back at least 
as far as ^Aiz93j . 

5. Expansion of the fundamental solution 

In this section we extend the results of Section 6 of [MNTOS] from the case 
71 = 3 to the case of arbitrary n > 3. 

We apply the method of bases with double orthogonality to the Cauchy problem 
(|4.1[) in the particular case, where A:' is a part of a ball B = B{0, R) with centre at 
the origin and radius R > 0. Let 5 be a smooth closed hypersurface in B which 
does not meet x — and divides B into two domains. Denote by X the closed 
domain that does not contain the origin. Its boundary dX consists of S and a part 
of the sphere dB in R", cf. Figure [TJ The advantage of using domains X of the 




Figure 1. A typical domain under consideration. 

above form lies in the fact that the problem reduces to analytic continuation from 
a small ball around to B. 

By the above, the double layer potential G{uo, fo) satisfies both A'PQ{uq, /q) = 
and {A + k^)0{uo, fo) = away from the closure of S. The latter equation is the 
scalar Helmholtz equation and follows from the former one. The idea is now to 
use bases with double orthogonality in Hilbert spaces of solutions to the Helmholtz 
equation to derive conditions of analytic continuation from B{Q, e) to B{0, R) for 
GiuoJo). 

In spherical coordinates, the Helmholtz operator A + k'^ in the space R" takes 
the form 

^ + k' = ^{{r^y + in -2)rl + - A,.-.) , (5.1) 

where /ign-i is the Laplace-Beltrami operator on the unit sphere. Recall that k is 
an arbitrary complex number with 3fc > 0. 

To solve the homogeneous equation {A + k'^)u = we use the Fourier method of 
separation of variables. Writing u{r,uj) — g{r, k)h{Lo) we get two separate equations 
for g and /i, namely 

(('^|:)' + ("-2)'^|: + ^V).9 ^ eg 

Agn-ift, — ch, 

c being an arbitrary constant. 

The second equation has nonzero solutions if and only if c is an eigenvalue of 
Z\s„-i. These are weh known to be c = j/(j/ + n - 2), for ly — 0,1, . . ., cf. |TS72| 
and elsewhere. The corresponding eigenfunctions of Ag^-i are spherical harmonics 
hy(ijj) of degree v, i.e. 

A^^-ihy ^ v{v + n-2)K. (5.2) 
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(( 



r— +{n-2)r— + {k'^r^-iy{iy + n-2)))g{r,k)^0. (5.3) 



Consider now the ordinary differential equation with respect to the variable 

r > 

d +("-2).- 

This is a version of the Bessel equation, and the space of its solutions is two- 
dimensional. 

For example, if n = 3 and k = 0, then g{r,0) = ar^ + hr^^^^ with arbitrary 
constants a and 6 is a general solution to (|5.3[) . In this situation any function 
r^hjj{Lp) is a homogeneous harmonic polynomial. In the general case the space of 
solutions to (j5.3p contains a one-dimensional subspace of functions bounded at the 
point r = 0, cf. |TS72| . 

For = 0, 1, . . ., fix a non-zero solution .g,y(r, k) of (j5.3p which is bounded at 
r = 0. Then 

(Z\ + fc2) (5.(r,fc)/i,.(^)) -0 (5.4) 

on all of M". Indeed, by ()5.1|) . (15. 2p and (15. 3|) we conclude that this equality holds 
in M" \ {0}. We now use the fact that gu{r, k)h^{uo) is bounded at the origin to see 
that llO)) holds. 

It is known that there are exactly J{i') linearly independent spherical harmonics 
of degree where 

^^''> ' in-2)M 

Pick an orthonormal basis 

i=i,. ..,./(..) 

in i2(§n-l). 

Lemma 5.1. For every R> 0, the system 

{by\r,u;,k):=g,{r,k)hi^\cj)} .^o,i,-.. (5.5) 

j=i,-.-, 

is an orthogonal basis in the subspace of L^{B{0, R)) consisting of solutions to the 
Helmholtz equation [A + fc^)w — 0. 

Proof. Indeed, as is an orthonormal basis in the space L^(S"~^) on 

the unit sphere, the system (|5.5p is orthogonal in L^{B{0, R)) because 

(6«,6(/))^.(s(o,fl)) = (/iW,4^))i2(s.-i) r g^{r,k)^jV;k)r^^~^dr 

Jo 

= 

for fi ^ V or i ^ j. Finally, since the system of harmonics {h\i^} is dense in 
C°°(S"^^) we see that system (|5.5p is dense in the subspace of L^{B{0,R)) con- 
sisting of the solutions of the Helmholtz equation in the ball, which completes the 
proof. □ 

For any fixed y e R" with \y\ > R, the fundamental solution e{x — y) of the 
Helmholtz equation satisfies {A + fc^)e(- — y) = in B{0, R) and is obviously square 
integrable in the ball. Therefore, e(- — y) can be represented in B{0, R) by its 
Fourier series with respect to the orthogonal system ()5.5p 

e{x -y)^Y.T. ^^'^(2^' ^) ^^'H^' (5.6) 
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the series converges in the L^{B{0, R)) -norm. The Fourier coefficients ci^\y,k) 
are defined by famihar formulas through the scalar product in L^(B{0, R)). These 
formulas show that c'J^y^k) are real analytic functions of y provided |?;| > R. 
Moreover, they satisfy the Helmholtz equation in the complement of B{0,R). On 
the other hand, the system (|5.5p is orthogonal in any space L^{B{0, R)) with R > 0. 
Hence, the coefficients c^^ (y, k) are actually independent of the particular choice 
of R satisfying R < \y\. This shows that c"J\y,k) are solutions of the Helmholtz 
equation away from the origin in R". Since e{x — y) = 0{\x — as x ^ y, we 

may tend R \y\ in the formulas for cji\y,k). This immediately yields explicit 
equalities 

c(^)(y,fc) = (e(.-y),6(^)(.,fc)) / ^'\gAr.k)\\-''dr 

\ / L''(B(0.\v\)) I . n 



for v — 0,1, . . . and j ~ 1, . . . , J(j^). 



Lemma 5.2. In the cone {{x,y) e R" x R" : |x|/|y| < 1}, one has the Fourier 
series expansion 15. 6\} . where the series converges uniformly together with all its 
derivatives on compact subsets of the cone. 



Proof. The Fourier series expansion is a direct consequence of Lemma [5.11 
The uniform convergence on compact subsets of the cone presents a more delicate 
problem. It can be handled in a familiar way, for b^^ are solutions of the Helmholtz 
equation. □ 

6. Regularisation 

Substituting (|5.6p into the formula for the fundamental solution <^*(a; — y) of 
M* we obtain 

oo 

'P\x-y)=Y,'l>t{x,y) (6.1) 

i/=0 

where the series converges uniformly along with all its derivatives on compact sub- 
sets of the cone {{x, y) e R" x M" : |x|/|y| < 1}, and 

'i'Ax,y) - [ ^A'{d^) -{i/ik)A'{d^)A"{d^) ) Z^^"^ \y^'^)''u 

for V = 0,1, . . .. We obtain the same formulas if we substitute ~~dy for the deriva- 
tives dx. 

Lemma 6.1. Every term 'I>^{x,y) is a real analytic matrix-valued function on 
the set R" x (R" \ {0}), satisfying 

M-{dx)K{x,y) = 0, 
NP'{dy){<l>l{x,y)f = 0. 

Proof. These properties are obvious by the very construction of degenerate 
kernels <Pl{x,y). The singularity at ?/ = is due to 

'''|5.(r,fc)pr"-idr. 

□ 
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A series expansion like (|6.ip with terms 4>l satisfying the transposed equa- 
tion M' {dy){<!>1{x,y))^ = is aheady sufficient to derive an exphcit formula for 
solutions of the Cauchy problem (|4.ip . Set 

N 

for {x, y) G M" x (M" \ {0}). In this way we obtain what is referred to as a Carleman 
function of the Cauchy problem, cf. Tar95, 10.4]. 

Theorem 6.2. Suppose {u, f) is a solution to Maxwell's equations in the inte- 
rior of X continuous up to the closure of S. Then 

for allxeX\ dX. 

Proof. We can assume, by approximating X by smaller domains whose bound- 
aries intersect the boundary of X only in 5, that (m, /) is smooth in X . Combining 
Corollary 12.31 with Stokes' formula and Lemma WT^ we conclude easily that 

for each fixed x ^ X \ dX and all = 0, 1, . . .. Let N oo. Since the series ()6.ip 
converges on dB{0, R) uniformly in x on compact subsets of -8(0, R), it follows that 
the part of the boundary integral over dX \ S tends to zero. This establishes the 
desired formula. □ 

Let (m, /) be a solution of the Cauchy problem (|4.ip . A straightforward com- 
putation shows that 

/ ^-(^'-^f rryi ,)ds-giuojo){^)-('!;i'\) (6.2) 

Js V«c^(^)Hhw)/ \Fn{x)J 

for X ^ S, where 

Write e > for the distance between S and the origin. If x S 5(0, e) then the 
left-hand side of (j6.2|) tends to zero, for the series (|6.ip converges uniformly on S. 
It follows that the series (i.e. a sequence of partial sums) {{Um^Fm)} converges 
to Q{uq^ /o) uniformly together with its derivatives on compact subsets of the ball 
B(0,e). 

The last observation and Theorem 14.11 fit together to yield certain conditions 
of solvability for the Cauchy problem. 

Corollary 6.3. If the series {{Un, Fn)} converges uniformly on compact sub- 
sets of the ball B(0,R), then the Cauchy problem ^T7p is solvable. 

Proof. Since the terms of the series {{Un, F^)} are component-wise solu- 
tions of the Helmholtz equation, it follows by the Stieltijes-Vitali theorem that its 
sum {U,F) = lim(t/Ar, Fjv) satisfies component-wise the same equation in 5(0, i?). 
Hence, {U, F) is a real analytic section of © on B{0, R). As ([/, F) actually 
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coincides with G{uo,fo) in the smaller ball B{0,e), the solvability of the Cauchy 
problem follows from Theorem 14. 11 as desired. □ 

In many interesting cases the uniform convergence of the series {{Un, F]\j)} is 
not only sufficient but also necessary for the solvability of (|4.ip . cf. [Shl96[ 2.9]. 



7. Propagation of electromagnetic waves 

We now turn to classical Maxwell's equation in a three-dimensional space, this 
latter can be a three-dimensional manifold X' as well. If X' is compact and closed, 
the theory is especially instructive, for the so-called Bohr-Sommerfeld radiation 
conditions are no longer needed. To demonstrate our constructions along more 
classical lines, we consider the case X' = M"^. As mentioned in Section [U classical 
Maxwell's equations have the form 

ikE + d*H = 0, 
-ikH + dE = 0, 

E and H being functions in a closed domain A" C M"^ with values in R'^. HE is 
suitably specified within 1 -forms and H within 2 -forms, both the exterior derivative 
d and its formal adjoint d* can be identified with the operator curl on vector fields 



Applying Corollary [273] to classical Maxwell's equations we obtain the Stratton- 
Chu formula [Str41) . 

Theorem 7.1. Suppose {E,H) is an electromagnetic wave in X whose electric 
component E and magnetic component H are both continuous up to the boundary. 
Then 

{l/ik)d*d -d* \ f -1 exp{ik\x ~ -l) f in{H) \ _ ( E{x) 



-d -{l/ik)dd* J Jg^ in \x--\ \-vM{E)) \H{x 
for all X G X \ dX , and the left-hand side vanishes away from X . 

Let X he a domain in M.^^ and S an open piece on dX of the form displayed in 
Figure [TJ We consider the problem of finding an electric field E and a magnetic 
field H in X with given tangential component Eq of E and normal component Hq 
of H on iS. 

Since the fundamental solution of the Helmholtz equation in M.^ is square inte- 
grable, the Fourier series expansion (|5.6p is evident. The coefficients c\^\y,k) are 
explicitly given by 

c\^\y,k)^(^ ^''''f^-~y'^\ b\]\;k)) J t'\g^{rM'r'dr 

\An \-~y\ J L^(B(o,\v\)) I 7o 

for i/ = 0, 1, . . . and j = 1, . . . , J{v), where J{v) = 2v + 1. 

Now the theory of Section [6] gives a Carleman function R]s[{x,y) of the Cauchy 
problem in explicit form. Namely, 

u=0 j = l 

the differential operator on the right-hand side acting in the variable x. 
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Theorem 7.2. Let {E,H) be an electromagnetic wave in X continuous up to 
S. Then 

\H{x)) N^^Js \-vM{E)) 

for all xeX\dX. 

In [Tar95| . several approaches to the Cauchy problem for solutions of linear 
elliptic equations with data on a part of the boundary are elaborated. They give 
not only explicit formulas of Carleman type for solutions but also conditions on 
the Cauchy data which are necessary and sufficient for the Cauchy problem to be 
solvable. The results demonstrate rather strikingly that the Cauchy problem for 
elliptic equations is overdetermined. It is solvable for a thin set of Cauchy data, 
and so any error in the Cauchy data leads to unsolvability. For this reason the 
variational approach to the Cauchy problem for elliptic equations is of practical 
interest. It works also in the case of nonlinear elliptic equations, cf. [LT09| . More 
generally, the problem consists in introducing reasonable classes of approximate 
solutions. In this sense any Carleman formula like the formula of Theorem 17.21 
produces approximate solutions to the Cauchy problem with data on a boundary 
piece for Maxwell's equations. 
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